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ABSTRACT 

The solar wind magnetic field contains rotations at a broad range of scales, which have been 
extensively studied in the MHD range. Here we present an extension of this analysis to the 
range between ion and electron kinetic scales. The distribution of rotation angles was found to 
be approximately lognormal, shifting to smaller angles at smaller scales almost self-similarly, 
but with small, statistically significant changes of shape. The fraction of energy in fluctuations 
with angles larger than a was found to drop approximately exponentially with a, with e- 
folding angle 9.8° at ion scales and 0.66° at electron scales, showing that large angles (a > 
30°) do not contain a significant amount of energy at kinetic scales. Implications for kinetic 
turbulence theory and the dissipation of solar wind turbulence are discussed. 
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1 INTRODUCTION 

The solar wind magnetic field displays broadband fluctuations that 
are in many ways consistent with our current understanding of 
plasma turbulence (e.g., Horbury, Wicks & Chen 2012; Bruno & 
Carbone 2013; Alexandrova et al. 2013). It also contains many 
small-scale structures, which may or may not arise from the tur¬ 
bulent dynamics. Determining the properties of these kinetic scale 
fluctuations is important for understanding the turbulent cascade 
and how its energy is dissipated, as well as the general structure of 
the solar wind. 

Large changes in the solar wind magnetic field direction are 
sometimes called “directional discontinuities” (Burlaga 1969) and 
are frequently (e.g., Burlaga 1971; Neugebauer et al. 1984; Tsu- 
rutani et al. 1996; Horbury et al. 2001; Paschmann et al. 2013) 
classified as rotational discontinuities (RDs) or tangential disconti¬ 
nuities (TDs) depending on which of these discontinuous solutions 
of ideal magnetohydrodynamics (MHD) (Landau & Lifshitz 1960) 
they more resemble. Their origin is debated (Burlaga 1969, 1971; 
Bruno et al. 2001; Vasquez et al. 2007; Borovsky 2008; Greco et al. 
2008, 2009; Owens, Wicks & Horbury 2011; Zhdankin et al. 2012; 
Zhdankin, Boldyrev & Mason 2012; Borovsky 2012; Malaspina & 
Gosling 2012; Arnold et al. 2013), in particular, whether they are 
generated in situ or represent plasma boundaries arising from pro¬ 
cesses at the Sun. For example, it has been proposed (e.g., Bruno 
et al. 2001; Borovsky 2008) that large angle changes represent flux 
tube boundaries in a filamentary picture of the solar wind originat¬ 
ing from early solar energetic particle observations (Bartley et al. 
1966; McCracken & Ness 1966). However, it has also been sug¬ 
gested that their waiting times and angular distributions are con¬ 
sistent with MHD turbulence (Vasquez et al. 2007; Greco et al. 


* E-mail: christopher.chen@imperial.ac.uk 


2008, 2009; Zhdankin et al. 2012; Zhdankin, Boldyrev & Mason 

2012) , suggesting that the majority are generated in situ. A recent 
simulation found MHD turbulence to produce both RDs and TDs, 
although RDs were found to be more numerous (Zhang et al. 2015). 

It is perhaps more instructive to examine the full distribution 
of magnetic field rotation angles, rather than just the large ones. 
Borovsky (2008) fitted the full distribution to a double exponential, 
interpreting the one dominating at small angles to be due to turbu¬ 
lence and the one at large angles to be due to flux tube boundaries. 
Borovsky (2012), however, later showed that plasma boundaries 
contribute only a small fraction of the distribution at all angles. Zh¬ 
dankin et al. (2012) showed that rather than an exponential, MHD 
turbulence simulations produce an angle distribution similar to the 
full distribution in the solar wind. Zhdankin, Boldyrev & Mason 
(2012) then showed that the solar wind distribution could be fit to 
a single lognormal population, and moreover that the shape of this 
is independent of scale over the three decades of the MHD inertial 
range, suggesting an underlying universal description. 

One of the questions addressed in this Letter is how the rota¬ 
tion angle distribution continues to below ion kinetic scales. Here 
the spectra of magnetic and density fluctuations steepen (e.g., Cole¬ 
man 1968; Russell 1972; Leamon et al. 1998; Smith et al. 2006; 
Chen et al. 2013b; Safrankova et al. 2013; Bruno & Trenchi 2014; 
Bruno, Trenchi & Telloni 2014; Chen et al. 2014a) and form a 
power law range with spectral index ~ —2.8 (Denskat, Beinroth 
& Neubauer 1983; Kiyani et al. 2009; Chen et al. 2010a; Alexan¬ 
drova et al. 2012; Chen et al. 2012; Sahraoui et al. 2013; Chen 
et al. 2013a). While this is generally consistent with scaling theo¬ 
ries of strong sub-ion scale turbulence (Vainshtein 1973; Biskamp, 
Schwarz & Drake 1996; Cho & Lazarian 2004; Schekochihin et al. 
2009; Chen et al. 2010b; Boldyrev & Perez 2012; Boldyrev et al. 

2013) , this range is less well understood than the turbulence at 
MHD scales. At kinetic scales the shapes of the distributions of 
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magnetic field component and density fluctuations do not vary 
much with scale (Kiyani et al. 2009; Wu et al. 2013; Chen et al. 
2014b), but angle distributions have not previously been measured. 

A related question is what type of structures the large an¬ 
gle changes correspond to. MHD turbulence is thought to generate 
sheet like structures, and both statistical (e.g., Biskamp 2003) and 
dynamical (Howes 2015) models have proposed to explain this. Re¬ 
cent work has shown that kinetic scale turbulence may do the same, 
for example, 3D fluid (Boldyrev & Perez 2012) and kinetic (Ten- 
Barge & Howes 2013) simulations show the development of 2D 
sheets at electron scales, although ID filaments have also been re¬ 
ported (Meyrand & Galtier 2013). It has been suggested that some 
of the current sheets formed in MHD turbulence (e.g., Carbone, 
Veltri & Mangeney 1990; Cowley, Longcope & Sudan 1997; Ser- 
vidio et al. 2009) and sub-ion scale turbulence (Haynes, Burgess & 
Camporeale 2014) may be undergoing reconnection. Evidence of 
small scale reconnection has been reported in the magnetosheath 
(Retino et al. 2007; Sundkvist et al. 2007) and solar wind (Gosling 
2012; Perri et al. 2012; Osman et al. 2014; Xu et al. 2015), possibly 
as part of the turbulent cascade. It has been suggested that a signif¬ 
icant amount (Leamon et al. 2000) or even the majority (Matthaeus 
et al. 2003) of energy in MHD turbulence is dissipated via recon¬ 
nection. However, others have suggested that enhanced dissipation 
can occur at current sheets through other processes such as Landau 
damping (TenBarge & Howes 2013). 

In this Letter, we extend the analysis of the distribution of 
magnetic field rotations into the kinetic range between ion and elec¬ 
tron scales. The scale-dependence of the distribution is examined, 
as well as the fraction of energy contained at the different angles. 
The implications for dissipation of turbulence and reconnection are 
also discussed. 


2 DATA SET 

The analysis was performed on a 7 hour interval of data from the 
Cluster spacecraft (Escoubet, Fehringer & Goldstein 2001) starting 
at 2003-02-12 23:00:00, when the four spacecraft were in the free 
solar wind (out of Earth’s foreshock). The typical spacecraft-frame 
correlation time at 1 AU is ~ 30 mins (e.g., Osman et al. 2014) 
so the interval covers many turbulence correlation lengths. For the 
magnetic field B, data from the DC magnetometers (FGM) (Balogh 
et al. 2001) and AC magnetometers (STAFF) (Cornilleau-Wehrlin 
et al. 2003) were combined into a single data set of 0.04 s resolution 
for each spacecraft, following the procedure given in Chen et al. 
(2010a). The spectrum of magnetic fluctuations remains at least 
an order of magnitude above the nominal instrumental noise floors 
over the entire range. 

To determine the plasma microscales, additional particle mo¬ 
ment data is required. The electron density n e was obtained from 
the plasma frequency determined from the high frequency electric 
field spectrum analyzer (WHISPER) (Decreau et al. 2001) and the 
perpendicular electron temperature Tj_ e from the electron electro¬ 
static analysers (PEACE) (Johnstone et al. 1997). Since data from 
the appropriate ion electrostatic analysers were not available, the 
ion (proton) velocity v; and perpendicular temperature 7j_ i were 
obtained from the Faraday cups (SWE) (Ogilvie et al. 1995) on 
the upstream Wind spacecraft (Acuna et al. 1995), shifted in time 
so that magnetic field features measured by both spacecraft were 
aligned. The mean parameters for this interval are: |B| = 10 nT, 
n e = 12 cnL 3 , |v;| = 390 km s -1 , = 5.6 eV, 7j_ e = 13 eV; the 

total perpendicular plasma beta is = 0.88. 




Figure 1. (a) PDFs of magnetic field rotation angle a for different time lags 
T. (b)PDFs of a = a/ (a), (c) Lognormal fits of PDFs of a; the inset shows 
the measured PDF (blue) and fit (red) for T = 0.04 s. The distribution for a 
Gaussian process is also shown (grey dashed). 


3 RESULTS 


The magnetic field rotation angle over time scale T and at time t is 


a(t, t) = cos 


B(f) -B(f+ t) 
|B(OI|B(f + T)| 


(1) 


Fig. la shows the probability density functions (PDFs) of a at 
scales t from 0.04 s to 82 s. Data from spacecraft 2^1 are included 
in the PDFs (spacecraft 1 data have significant spin tone). The dis¬ 
tribution moves to smaller angles at smaller scales, as expected, 
since the fluctuation amplitudes decrease. 

Under the Taylor (1938) hypothesis, spacecraft-frame time 
scales correspond to spatial scales in the plasma. While this may 
break down if the turbulence evolves quickly as it passes the space¬ 
craft, it is thought to be a good approximation for anisotropic ki¬ 
netic Alfven turbulence (Howes, Klein & TenBarge 2014; Klein, 
Howes & TenBarge 2014), which measurements suggest is the 
dominant component at kinetic scales (Chen et al. 2010a, 2013a; 
Telloni, Bruno & Trenchi 2015). Under this approximation, t ~ 
1.3 s corresponds to ion scales (e.g., gyroradius, inertial length) 
and T ~ 0.04 s to electron scales. The data in Fig. la show that the 
maximum value of a in this interval is 65° at ion scales and 5.4° at 
electron scales; the mean values are 2.2° and 0.17° respectively. 

To test the self-similarity of a over the kinetic range, Fig. lb 
shows the distributions of a = a/ (a), i.e., a normalised to its 
mean value, where the angular brackets denote an average over t. 
Apart from variations at large angles for large T, the distributions 
appear to be quite self-similar (although non-Gaussian) with only 
small changes with scale. However, a quantitative test is required to 
assess the degree of self-similarity. The kurtosis of the a distribu¬ 
tion was measured directly, although this was found to be unreliable 
due to the finite interval length (Horbury & Balogh 1997; Dudok de 
Wit 2004). Therefore, the distributions were fit to various functions. 
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Figure 3. PDFs of % (Equation 3) at different time lags T. Colors are the 
same as those in Fig. 1. 


Figure 2. Lognormal fit parameter fi and other inferred parameters of the 
distribution of a as a function of time lag T. Scales corresponding to the pro¬ 
ton (dashed) and electron (dash-dotted) gyroradii (red) and inertial lengths 
(green) are marked, along with the measured spectral break (black dashed). 


The best fit was obtained for a lognormal distribution, 


/(*) = 


i 


\a\f2K 


exp 


(liuc — /l ) 2 


2c 2 


( 2 ) 


similar to previous findings for the MHD range (e.g., Bruno et al. 
2004; Zhdankin, Boldyrev & Mason 2012). The fitted distributions 
are shown in Fig. lc; they are able to capture most of the features 
of the measurements, although the smallest angles are slightly un¬ 
derestimated (see inset). The fit is better than a double exponential 
(used by Borovsky 2008), and has fewer free parameters. 

Fig. 2 shows the fit results as a function of T. Since at all 
scales the mean of a is 1, the fit effectively has one free parameter, 
so only fi is shown. The mode, skewness and kurtosis calculated 
from the fit parameters are also shown. The error bars represent 
the standard deviations calculated from 7 non-overlapping 1-hour 
sub-intervals. It can be seen that there is a small but statistically 
significant variation in the parameters with scale. For example, the 
kurtosis increases from large to small scales, peaks at ion scales, 
then decreases towards electron scales; a behaviour consistent with 
the kurtosis of the magnetic field component fluctuations in this 
interval (not shown here) and other works (e.g., Wu et al. 2013). 

Zhdankin. Boldyrev & Mason (2012) used the parameter 


XM- ' SB/B ~ s l% W2) ' . 0 ) 

where 8B/B = |B(f + t) — B(f)|/|B(?)|, to measure the degree to 
which the magnetic field changes magnitude during a rotation. 
X = 0 corresponds to a pure rotation and 0 < % ^ 1 corresponds to 
a rotation with change in field strength (note that this does not dis¬ 
tinguish TDs from RDs). Fig. 3 shows the PDFs of X over the same 
range of scales as Figures 1 and 2. At all scales the distribution 
peaks at x = 0 showing that the fluctuations are predominantly pure 
rotations. However, the peak becomes less strong towards smaller 
scales, consistent with the turbulence being more compressive in 
the kinetic range. Since the fluctuations are mostly pure rotations, 
8B/B ps a for small a , and indeed the results for the distributions 
of 8B/B match the above results for the distributions of a. 

Finally, to help understand the distribution of energy and how 
it is dissipated at kinetic scales, the fluctuation energy contained in 
the magnetic rotations was measured. Fig. 4 shows the fraction of 


angles larger than a, and the fraction of magnetic fluctuation en¬ 
ergy |B(f + t) — B(f)| 2 that they contain, at ion and electron scales. 
For both scales, the energy drops approximately exponentially with 
a. A similar result was found for the energy dissipation in current 
structures in a recent 3D RMHD simulation (Zhdankin et al. 2014), 
although the comparison is less good with a 2D PIC simulation 
(Wan et al. 2012), perhaps suggesting that 3D simulations are bet¬ 
ter able to capture the properties of turbulent structures formed in 
the solar wind. The e-folding value is 9.8° for ion scales and 0.66° 
for electron scales, showing that throughout the kinetic range there 
is only a small amount of energy in large angles (a > 30°), al¬ 
though since the typical angles are small, those a few times larger 
than the mean still contain a significant fraction of the energy. 

The above analysis was performed on several other intervals 
of data and found to be consistent with that presented here. 


4 DISCUSSION 

We have shown that the distribution of magnetic field rotations is 
approximately lognormal from ion to electron scales, however there 
is a small, but statistically significant, change in shape with scale, 
e.g., the kurtosis decreases by a factor of 1.6 from ion to elec¬ 
tron scales, consistent with the magnetic field components (e.g., 
Wu et al. 2013). The interval is not long enough to probe the full 
MHD range to enable a full comparison with the results of Zh¬ 
dankin, Boldyrev & Mason (2012), although the lognormal fit pa¬ 
rameters are somewhat different in the range where the analyses 
overlap (3 s to 82 s). This may be related to the fact that the cur¬ 
rent study is based on 7 hours of data (a single “stream”), while 
the study of Zhdankin, Boldyrev & Mason (2012) used 6 years, 
meaning that the distributions there may in part reflect the varying 
driving conditions rather than being solely an effect of the turbulent 
cascade. 

Figures 1 and 4 also show that the rotation angles at kinetic 
scales are relatively small, much smaller than those obtained by 
Perri et al. (2012). The reason for this is that Perri et al. (2012) 
used only data from STAFF in their calculation of a (rather than 
the combined FGM-STAFF data set used here), which effectively 
provides a magnetic field time series high pass filtered at ~ 1 Hz. 
Due to the lack of DC field, the values of a obtained by Perri et al. 
(2012) are much larger than the true rotation angles, as shown here. 
The larger fluctuations of the distribution may still correspond to 
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Figure 4. (a) Fraction F of ion scale fluctuations with an angle > a (blue) 
and fraction of magnetic fluctuation energy in those fluctuations (red). The 
black dashed line shows an exponential fit F = 0.96exp(—a/9.8), (b) Same 
for electron scales with fit F = 0.73exp(—a/0.66). 

particular types of structures, although current measurements are 
not sufficient to unambiguously determine their geometry. 

Given the small amount of energy in large angle fluctuations 
at kinetic scales (Fig. 4), it might be questioned how significant re¬ 
connection is for dissipating turbulence in the solar wind. While 
this cannot be answered with the current data set, knowing the dis¬ 
tribution of angles can help constrain this possibility. For example, 
reconnection, in the traditional sense, is known to be suppressed 
at small angles if /3 fluctuations across the current sheet are large 
enough (Phan et al. 2010) due to the diamagnetic drift of the x-line 
being faster than the outflow speed (Swisdak et al. 2003, 2010). 
Fligh resolution, low noise particle data are required to check this 
condition for turbulent fluctuations at kinetic scales. 

Reconnection has been identified in the solar wind from 
Alfvenic jets associated with bifurcated magnetic field rotations 
(e.g., Gosling 2012), although typically only a few events per day 
are found overall (e.g., Gosling et al. 2007; Osman et al. 2014), with 
fewer than one per day in the purer fast wind (Gosling 2007). While 
reconnection at small angles is occasionally detected (Gosling & 
Phan 2013), this is only for very low /3 plasma not typical of the 
solar wind at 1 AU. The typical separation between measured re¬ 
connection events is one to two orders of magnitude larger than the 
turbulence correlation length and such rarity led Gosling (2007) to 
question whether it can be a major way for the turbulence to be dis¬ 
sipated. Fligh resolution magnetic field data have shown that only a 
small number of reconnection events (identified from bifurcated ro¬ 
tations) are missed due to limited resolution particle data (Gosling 
& Szabo 2008). The data presented here were examined for bifur¬ 
cated magnetic rotations, although no clear signatures were found. 
It may be the case that the amount of turbulent energy dissipated 
via reconnection (rather than dissipated directly) depends on pa¬ 


rameters of the system, such as the length of the inertial range or 
plasma |3. However, even if some energy conversion occurs through 
reconnection, how it is finally dissipated still remains to be identi¬ 
fied (e.g, Loureiro, Schekochihin & Zocco 2013; Haynes. Burgess 
& Camporeale 2014; Drake & Swisdak 2014). 

Finally, knowing the full distributions of magnetic fluctuations 
at kinetic scales allows one of the important assumptions of sev¬ 
eral approaches to solar wind turbulence to be tested. It is some¬ 
times assumed, for example in reduced non-linear models (e.g., 
Howes et al. 2006; Schekochihin et al. 2009; Boldyrev et al. 2013; 
Kunz et al. 2015) and linear theory (e.g., Gary & Smith 2009), that 
5B/B <C 1. While this is well satisfied for the typical fluctuation (on 
average 8B/B = 0.040 at ion scales and 8B/B = 0.0034 at electron 
scales for this interval), it has been questioned to what extent this is 
valid for the larger amplitude fluctuations which are an important 
part of the cascade. In the current interval, at ion scales 1.1% of the 
fluctuations, containing 30% of the energy, have 8B/B > 0.2 and at 
electron scales all fluctuations satisfy 8B/B < 0.1. Therefore, one 
could conclude that the assumption becomes valid for the major¬ 
ity of fluctuations around ion scales, and is very well satisfied by 
electron scales. 
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